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ABSTRACT 

We continue to investigate which polynomials can possibly 
occur as factors in the denominators of rational solutions of 
a given partial linear difference equation. In an earlier arti- 
cle we had introduced the distinction between periodic and 
aperiodic factors in the denominator, and we gave an algo- 
rithm for predicting the aperiodic ones. Now we extend this 
technique towards the periodic case and present a refined 
algorithm which also finds most of the periodic factors. 
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1.1.2 [Computing Methodologies]: Symbolic and Alge- 
braic Manipulation — Algorithms 

General Terms 

Algorithms 
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1. INTRODUCTION 

The usual approach for finding rational solutions of lin- 
ear difference equations with polynomial coefficients is as 
follows. First one constructs a nonzero polynomial Q such 
that for any solution y = p/q of the given equation we must 
have q \ Q. Such a polynomial Q is called a denominator 
bound for the equation. Next, the denominator bound is 
used to transform the given equation into a new equation 
with the property that a polynomial P solves the new equa- 
tion if and only if the rational function y = P/Q solves the 
original equation. Thus the knowledge of a denominator 
bound reduces rational solving to polynomial solving. 
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The first algorithm for finding a denominator bound Q 
was given by Abramov in 1971 [2, 3, 4]. During the past 
fourty years, other algorithms were found [12, 15, 6, 8] and 
the technique was generalized to matrix equations [1, 5] as 
well as to equation over function fields [13, 7, 14]. Last 
year [11] we made a first step towards a denominator bound- 
ing algorithm for equations in several variables (PLDEs). 
We found that some factors of the denominator are easier 
to predict than others. We called a polynomial periodic if 
it has a nontrivial gcd with one of its shifts, and aperiodic 
otherwise. For example, the polynomial 2n — 3fc is periodic 
because shifting it twice in k and three times in n leaves 
it fixed. We say that it is periodic in direction (3,2). An 
example for an aperiodic polynomial is nk + 1. The main 
result of last year's paper was an algorithm for determining 
aperiodic denominator bounds for PLDEs, i.e., we can find 
Q such that whenever y = -2- solves the given equation and 
q is aperiodic, then q \ Q. 

The present paper is a continuation of this work. We 
now turn to periodic factors and study under which circum- 
stances a slightly adapted version of last year's algorithm can 
also predict periodic factors of the denominator. We pro- 
pose an algorithm which finds the periodic factors for almost 
all directions. Every equation has however some directions 
which our algorithm does not cover. But if, for instance, 
we have a system of two equations and apply our algorithm 
to each of them, then the two bounds can under favorable 
circumstances (which can be detected algorithmically) com- 
bined to a denominator bound which provably contains all 
the factors that can possibly occur in the denominator of 
any solution of the system. This was not possible before. So 
while until now we were just able to compute in all situa- 
tions some factors, we can now also find in some situations 
all factors. 

Despite this progress, we must confess that our results 
are still of a somewhat academic nature because denomina- 
tor bounds in which some factors are missing are not really 
enough for solving equations. And even when a full denomi- 
nator bound is known, it still remains to find the polynomial 
solutions of a PLDE, and nobody knows how to do this — the 
corresponding problem for differential equations is undecid- 
able. But in practice, we can heuristically choose a degree 
bound for finding polynomial solutions, and knowing parts of 
the possible denominators is certainly better than knowing 
nothing, and the more factors we know, the better. Apart 
from this, we find it interesting to see how far the classi- 
cal univariate techniques carry in the multivariate setting, 



and we would be curious to see new ideas leading towards 
algorithms which also find the factors that we still miss. 

2. PREPARATIONS 

Let K be a field of characteristic zero. We consider poly- 
nomials and rational functions in the r variables m , . . . , n r 
with coefficients in K. For each variable rii, let Ni denote 
the shift operator mapping rii to n, + 1 and leaving all other 
variables fixed, so that 

Niq(m,. . . ,n r ) 

= g(m, . . . ,7ii_i,n» + 1, n i+ i,. . . ,n r ) 

for every rational function q. Whenever it seems appropri- 
ate, we will use multiindex notation, writing for instance n 
instead of m , . . . , n r or N l for A?] 1 N l 2 2 ■ ■ ■ Np~ . 
We consider equations of the form 

J2 a * NS y = f C 1 ) 

ses 

where S C ZT is finite and nonempty, / G K[n] and a s G 
K[n] \ {0} (s G S) are given, and y is an unknown rational 
function. Our goal is to determine the polynomials p G K[n] 
which may possibly occur in the denominator of a solution y, 
or at least to find many factors of p. 

We recall the following definitions and results from our 
previous paper [11]. 

Definition 1. Let p, q, d € K[n]. 

1. The set Spread(p,q) := { i G 7L r : gcd(p,N l q) / 1} 
is called the spread of p and q. For short, we write 
Spread(p) := Spread(p,p). 

2. The number Disp fc (p, q) := max{|ifc| : (ii,...,i r ) G 
Spread(p, q)} is called the dispersion of p and q with 
respect to k G {1, . . . , r}. (We set maxyl := — oo if A 
is empty and max A := oo if A is unbounded.) 

3. The polynomial p is called aperiodic if Spread(p) is 
finite, and aperiodic otherwise. 

4- The polynomial d is called an aperiodic denominator 
bound for equation (1) i/d / and every solution y 
can be written as for some a,b,u G K[n] where u 
is periodic and b \ d. 

5. A point p G S C ZT C ]R r is called a corner point of S 
if there exists a vector v G R r such that (s — p) • v > 
for all s G S \ {p}. Such a vector v is then called an 
inner vector, and the affine hyperplane H := {x G TR r : 
(x — p) ■ v = 0} is called a border plane for S. 




Theorem 1. Let p, q G K[n], 

1. If p is irreducible, then Spread(p) is a submodule ofTT 
and p is aperiodic if and only i/Spread(p) = {0}. 



2. If p and q are irreducible, then there exists s£F such 
that s + Spread(p, q) is a submodule of HT . 

3. There is an algorithm for computing Spread(p, q). 

4- There is an algorithm for computing an aperiodic de- 
nominator bound for (1) given the support S and the 
coefficients a a (s G S). 

3. DENOMINATOR BOUNDS MODULO A 
PRESCRIBED MODULE 

Our goal in this section is to determine the factors whose 
spread is contained in some prescribed set W C 2 r . Under 
suitable assumptions about W such factors must pop up in 
the coefficients of the equation (cf. Lemma 2 below) and 
under stronger assumptions we can also give a bound on 
the dispersion between them (cf. Theorem 2 below). Using 
these two results we obtain a denominator bound relative 
to W (cf. Theorem 3 and Algorithm 1) below. In the next 
section, we then propose an algorithm which combines the 
denominator bounds with respect to several sets W. It turns 
out that by considering only finitely many sets W one can 
obtain a denominator bound with respect to infinitely many 
sets W. 

Definition 2. Let W C 7L r with G W . A polynomial 
d G K[n] \ {0} is called a denominator bound of (1) with 
respect to W if for every solution y — p/q G K(n) of (1) 
and every irreducible factor u of q with Spread (u) C W we 
have u \ d. 

Typically, W will be a submodule of 7L r or a union of such 
modules. The definition reduces to the notion of aperiodic 
denominator bound when W = {0}. In the other extreme, 
when W = ZT then d is a "complete" denominator bound: 
it contains all the factors, periodic or not, that can possibly 
occur in the denominator of a solution y of (1). In general, 
d predicts all aperiodic factors in the denominator of a solu- 
tion as well as the periodic factors whose spread is contained 
in W. 

Denominator bounds with respect to different submodules 
can be combined as follows. 

Lemma 1. Let Wi, . . . , W m be submodules of HI , and let 
di,...,d m be denominator bounds of (1) with respect to 
W\, . . . , W m , respectively. Then d := lcm(cii, . . . , d m ) is a 
denominator bound with respect to W := W\ U • • • U W m . 

Proof. Let u be an irreducible factor of the denominator 
of some solution of (1) and suppose that U := Spread(u) C 
W . It suffices to show that then U C Wk for some k, because 
then it follows that u \ dk | d, as desired. 

We show that if U contains some vector x ^ Wi, then 
Wi C W 2 U ■ ■ ■ U W m , hence U CW 2 U---U W m - Applying 
the argument repeatedly proves that U C Wk for some k. 

Let y G Wi. Since U is a submodule of 7L r , we have 
x + ay G U for all a G Z. By assumption U C WiU- ■ -UW m , 
so each such x + ay must belong to at least one module We 
(£ = l,...,m). It cannot belong to Wi though, because 
together with y G W\ this would imply x G Wi, which 
is not the case. Therefore: For every a G 7L there exists 
I G {2, . . . , m } such that x + ay G We . 

Since 7L is infinite and m is finite, there must be some index 
I G {2, . . . , m} for which there are two different ol\ , a 2 G 7L 



with x + a?iy £ We and x + a^y € We- Since We is also 
a submodule of Z r , it follows that (ai — a^)y £ Wi, and 
finally y 6 Wi C W 2 U • • • U W m , as claimed. □ 

The next result says that factors of denominators tend to 
leave traces in the coefficients of corner points of S. 

Lemma 2. Let W be a submodule of IT and let u with 
Spread(-u) (- W be an irreducible factor of the denominator 
of some solution y of (1). Let p £ S be a corner point of 
S with an inner vector v £ !R r orthogonal to W (meaning 
w ■ v = for all w £ W ). Then there exists i £ TIT such 
that N'u | a-p. 

Proof. If u' is another irreducible factor of the denom- 
inator of y and Spread (u,u') is nonempty, then we have 
Spread(u, u) C c + W for some c £ 717 . This follows from 
Theorem 1.2 and the assumption Spread(u) C W . For the 
full denominator d of y, we can thus find ci, . . . , c m £ 7L r 
with Spread(«, d) C [JStLi ( c fe + W) where each element from 
C is necessary. Let C = {ci , . . . , c m } be such a choice, and 
let i £ {1, ... , m} be such that Ci ■ v is minimal. 

We have 

a P N p y = f- a * NB y 
ses\{ P } 

as an identity in K(n). Therefore, every factor in the de- 
nominator of N p y must either be canceled by a p or it also 
occurs as a factor in at least one of the N a y (s £ S \ {p})- 
The factor N p+Ci u appears in the denominator of N p y. 
If it also appeared in the denominator of N s y for some 
s £ S \ {p}, then this would imply N p+C 'u = N s+C ^u for 
some j £ {1, . . . , m}. But then s — p + Cj — a £ W, which 
is in contradiction to 

(s - p + Cj - Ci) ■ V = (s - p) ■ V + (cj - Ci) ■ V / 

>0 >0 

because W is orthogonal to v by assumption. Hence N p+Ci u 
cannot appear as a denominator on the right hand side, and 
hence it must be canceled on the left hand side. This forces 
7V P+Ci it | <z p , so the claim is proven for i := p + Cj. □ 

The lemma tells us for which choices of W C 7L r some- 
thing nontrivial may happen. Let us illustrate this with an 
example. 

Example 1. The equation 
(Ik -2n + l)(k + n + l)y(n, k) 

+ (8k 2 + 2kn + k + 6n 2 + 13n + 6)y(n, k + 1) 

- 2(6k 2 + 2kn + 13k + 2n + n + 6)y(n + 1, k) = 

has the solution y = (n 2 + 2k 2 ) /(k + n + 1). Its denom- 
inator is periodic, Spread(fc + n + 1) = (J^ZI. Lemma 2 
predicts the appearance of k + n + 1 (or at least some shifted 
version of it) in the coefficient of y(n,k), because for the 
choice W = (^-JZj, the point p = (?) £ S admits the choice 
v = (V) in accordance with the requirements imposed by the 
lemma. Note that no shift equivalent copy of k + n+1 occurs 
in the coefficients of y(n,k + 1) or y(n + l,k), which does 
not contradict the lemma, because the points (A and (J) lie 
on a line parallel to W . This has the consequence that for 
these points, there does not exist a vector v with the required 
property. 



Conversely, the factor 4fc — 2n + 1 cannot possibly ap- 
pear in the denominator of a solution, because for W' := 
Spread(4fc — 2n + 1) = ( 2 )Z we can take p' = (J) and 
v' = (^y 2 )> an d according to the lemma, some shifted ver- 
sion of 4fc — 2n + 1 would have to appear in the coefficient 
ofy(n + l,k). 

More generally, for any nontrivial submodule W" of 
other than W , Lemma 2 excludes the possibility of periodic 
factors whose spread is contained in W" , because such fac- 
tors would have to leave a trace in at least one of the coeffi- 
cients of the equation. 




In the previous example, we could thus determine all the 
interesting modules W by just looking at the spreads of the 
the factors of the coefficients of the equation. The following 
example indicates that this is not always sufficient. 

Example 2. The equation 

(2k - 3n 2 - 8n - 5)y(n, k + 1) 
+ (k + 3n 2 + 5n + 4)y(n + 1, k) 

- (5k - 3n 2 - lln - 7)y(n + 1, k + 1) 
+ (2k - 3n - 8n - 3)y(n + 2, k) = 

also has the solution y = (n 2 + 2k 2 )/ (k + n+1). Its denom- 
inator k + n + 1 does not appear in any of the coefficients of 
the equation. This is because for its spread W = there 
are no suitable p £ S and v £ K 2 matching the conditions 
of the lemma because the points (J), (Tj as well as the points 

(o)' (i) ^ e on a ^ ne P ara ^d t° S. 

In summary, in order for W to be the spread of a factor 
that can appear in the denominator of a solution of (1), 
W must be contained in the spread of some coefficient of the 
equation (as in Ex. 1) or it must be parallel to one of the 
faces in the convex hull of the support S (as in Ex. 2). For 
every equation, we can thus determine some finitely many 
submodules of of codimension one such that each possibly 
occuring spread W is contained in at least one of them. 

3.1 A Normalizing Change of Variables 

Let 7L r = V © W be a decomposition of TIT into sub- 
modules. Our goal is to obtain denominator bounds with 
respect to W by applying the algorithm from last year [11] 
to V = T7 /W. It turns out that this can be done provided 
that W is sufficiently nondegenerate. In order to formulate 
the precise conditions on W without too much notational 
overhead, it seems convenient to make a change of coordi- 
nates. 

Let invertible matrices A = ((oi,j))i J=1 £ (Q rxr act on 



K(n) via 

A ■ 2/(m, • • • ,n r ) := y(a 1A n 1 + a lj2 n 2 + h ax tT n r , 

02,1^1 + a 2t2 n 2 + ■ ■ ■ + a2, r n r , 



ar,Wl + flr,2^2 + + a r ,r^r)- 

We obviously have A ■ (p + g) = (A • p) -f ( A ■ g) and A ■ (pg) = 
(A ■ p)(A ■ g) for all p, g G K(n). It can be checked that we 
also have 

A-(N s y) = N A " s (A-y) 

for every A G % TXr with | det = 1 and every s G Z r and 
every y G K(n). It follows that y G K(n) is a solution of 
(1) if and only if y = A -1 • y is a solution of the transformed 
equation 

J2(A-a s )N A " s y = A-f, 

ses 

or equivalently of 

J2&sN s y = f, 

ses 

where S = {A~ 1 s : s G S}, a B := A ■ a,As (s G S), and 
/ = ^ • /• 

Now take A G Z r with | det A| = 1 such that the first t 
rows of A form a basis of V and the last r — t rows of A 
form a basis of W . Then the transformation just described 
maps the basis vectors of V to the first t unit vectors and 
the basis vectors of W to the last r — t unit vectors. In 
other words, we can assume without loss of generality that 
V itself is generated by the first t unit vectors and W by the 
last r — t unit vectors in TIT . We will make this assumption 
from now on, unless otherwise stated. Note that this con- 
vention implies for an irreducible polynomial u G K[n] that 
Spread(u) = W is equivalent to u being free of the variables 
nt+i,n t +2, • ■ • , n r and aperiodic as element of K[ni, . . . ,n t ]. 

By applying, if necessary, a suitable power of N\ on both 
sides of the equation we can further assume without loss of 
generality that min{si : (si, . . . , s r ) G S} = 0, and we set 
k := max{si : (si, . . . , s r ) G S}. 

3.2 Bounding the Dispersion 

With this transformation w.r.t. the submodule W of 7L r 
and under the assumption that the extreme points (2) of 
S have certain properties, Theorem 2 explains how one can 
bound the dispersion along the x i -coordinate of all factors 
/ with Spread(/) G W that occur in the denominator of 
a solution of (1). This result is a refinement of Lemma 2 
from [11]. 

Lemma 3. Let it, v G K[n]\{0} with Spread(w) C W and 
Spread(v) C W. Then 

|{(s!,...,s t ) G %* I 3(s t+1 ,...,s r ) G Z'-' : 

N ( - Sl '-' 3r) u = v}\ < 1. 

Proof. Take s, s' with N s u = v = N s> u. As N s ~ s ' u = u, 
it follows s - s' G W = {0}' x Z r ~ t , and thus the first t 
components of s, s agree. □ 



Theorem 2. Let 

A = {( Sl ,...,s r ) G5 :S i=0}, 

B = {{si,...,s r )€S:s 1 = k}. ( ' 

Suppose that no two elements of A agree in the first t co- 
ordinates, and that the same is true for B. Let a[ be those 
polynomials which contain all irreducible factors f of ai with 
Spread(/) C W. Let 

s := max{Disp 1 (ag, N± k a' t ) : s G A and t G B}. 

Then for any solution y = p/q G K(n) of (1) and any irre- 
ducible factors u,v of q with Spread(u), Spread(w) C W we 
have Disp 1 (ti, v) < s. 

Proof. As S is not empty, A,B axe nonempty. W.l.o.g. 
we may assume that the minimal element of A w.r.t. lexico- 
graphic order is the zero vector. 

Suppose that there are irreducible factors u,v of q with 
Spread(w) C W, Spread(u) C W and d := Disp 1 (it,i)) such 
that d > s; take such u, v such that d is maximal. Consider 
all the factors N u u and N v v occurring in q where the first 
entry in u and v is 0. Note that by Lemma 3 there are only 
finitely many choices of the first t components, so we can 
choose two such factors from q where the first t components 
of u are minimal and the first t components of v are maxi- 
mal w.r.t. lexicographic order; these factors are denoted by 
u',v' respectively. 

• First suppose that v! divides one of the polynomials a s 
with s G A. In this case we choose the polynomial a w with 
w = (wi, . . . , w r ) G B such that (w 2 , . . . ,wt) is maximal 
w.r.t. lexicographic order (uniqueness is guaranteed by the 
assumption that no two elements from B agree in the first t 
components). We can write (1) in the form 

ATy=i(/- J2 a * NS y)- ( 3 ) 
w ses\{w} 

Now observe that the factor N™v' does not occur in the 
denominator of any N s y with s G S \ {w}: 

1. Suppose that there is s G S\B such that iV w i/ occurs 
in N B q, i.e., N™~ s v' is a factor of q. Since the first 
component of w is k (w G B) and the first component 
of s is smaller than k (s ^ B), the first component of 
w — s is positive. Moreover, since the distance between 
the factors v' and u' of q is d in the first component, 
the factors v' and iV w ~ s i/ of q have distance larger 
than d in the first component; a contradiction that the 
distance d is maximally chosen. Consequently, if A rw w' 
is a factor in the denominator of N s y with s € S, it 
follows that s G B. 

2. Suppose that there is s G B with w ^ s such that 
N™v' is a factor of N s q. Then A W ~V is a factor of 
q. Since the first component of the vectors in B is 
k, but the first t components in total cannot be the 
same for two different vectors of B, it follows that the 
first entry in w — s is zero and at least one of the 
others is non-zero; in particular, by the maximality 
assumption on w the first non-zero entry is positive. 
Hence we find v' = (0,v' 2 , ■ ■ ■ , v' r ) :— v + w — s such 
that N v v is a factor of q and such that (w 2 , . . . , v' t ) is 
larger than (v 2 , . . . ,v t ) w.r.t. lexicographic ordering; a 
contradiction to the choice of the vector v. 



Since /, a s G K[n], the common denominator of the rational 
function / — X] s es\{w} a aN s y does not contain the factor 
N™v'. Now suppose that N™v' is a factor of a w . Since 
w G B, its first component is k. But then, since u' and 
v' have distance d in the first coordinate, also the factors v! 
and N^ k N m v' have distance d. Thus Dispi (a s , jVf fc a w ) > d 
which implies that s > d; a contradiction. Overall, the com- 
mon denominator on the right hand side of (3) cannot con- 
tain the factor AT W ?/ which implies that the denominator of 
iV w y is not divisible by N™v' . Thus the denominator of y, 
in particular q is not divisible by v'; a contradiction. 
• Conversely, suppose that u' does not divide any of the 
polynomials a s with s £ A. Now let w = (0, 102 • • ■ , w T ) € A 
such that (u;2, . . . ,w t ) is minimal w.r.t. lexicographic order- 
ing (again it is uniquely determined by the assumptions on 
^4), and write (1) in the form (3); by our assumption stated 
in the beginning, w is just the zero vector 0. By analogous 
arguments as above (the roles of A and B are exchanged) 
it follows that u' does not occur in the denominator of any 
N s y with s G S \ {0}. Hence as above, the common denom- 
inator of / — X] s es\{o} a sN s y does not contain the factor 
u . Moreover, since u does not divide any a s from s G A, 
the factor u does not occur in do. In total, the factor it' is 
not part of the denominator on the right hand side of (3), 
but it is a factor of the denominator on the left hand side; 
a contradiction. □ 

If the required properties on the sets (2) in Theorem 2 are 
violated, our bounding strategy does not work, as can be 
seen by the following example. 

Example 3. Fix W := Spread(fc + n + 1) = (\)% and 
take V — The problem from Example 1 is normalized 

by the change of variables n — > k and k — > n — k (i.e., a basis 
transformation fj }^with determinant —1 is chosen) and 

one obtains V' = ( )Z and W' = This gives the new 

equation 

(n + l)(-6fe + An + l)y(n, k) 

+ (12k 2 - Unk + 12k + 8n 2 + n + 6)2/(71 + 1, fc) 
- 2(6k 2 - lOnfc - 12k + 6n 2 + 13n + 6)y(n + 1, k + 1) = 

with the new structure set S' = {( )j (q)i (1)} which now 
has the solution y = 3fc ~ 4nk + 2n where the denominator 
consists of the factor n+l with Spread(n + 1) = W' . As 
observed already in Example 1 one can predict the factor n+l 
(up to a shift in n) by exploiting Lemma 2. However, one 
cannot apply Theorem 2. For S' we get the sets A = { (°) } 
and B = {(J), (})} where in B the two vectors are the same 
in the first component but differ in the second component. 

3.3 Denominator Bounding Theorem 

The denominator bounding theorem says that if we rewrite 
the equation (1) into a new equation whose support contains 
some point p which is sufficiently far away from all the other 
points in the support, then we can read off a denominator 
bound from this new equation. We will need the following 
fact, which appears literally as Theorem 3 in [11] (with W, S' 
renamed to R~ , R + here in order to avoid a name clash with 
the meaning of W in the present paper). 



Lemma 4. Let p be a corner point of S with border plane 
H and inner vector v. Then for every s > there exist finite 
sets 

fl" C2 r n I) (H + ev) and 

0<e<s 

R + C Hi n [J(H + ev), 

and polynomials b,bi £ K[n] such that for any solution y £ 
K(n) of (1) we have 



N p y 



_ b + J2 ieR+ b^y 



UieR- Ni ~ Pa e 



(4) 



The sets R~ and R + and the polynomials b, bi can be com- 
puted for a given s, S, p, and v by Algorithm 2 from [11]. 
The next theorem provides a denominator bound with re- 
spect to W . It is an adaption of Theorem 4 from [11] to the 
present situation. We continue to assume the normalization 
V = x {0} r -*, W = {0}' x TL r -*. 

Theorem 3. Let s 6 NU{- 00} be such that for any solu- 
tion y = p/q G K(n) of (1) and any irreducible factors u, v 
of q with Spread(it), Spread(i)) C W we have Disp 1 (u, v) < 
s. Let p be a corner point of S for which there is an inner 
vector v = (yi, . . . , v r ) with vi > 1 as well as an inner vec- 
tor v' orthogonal to W. For these choices of s, p, and v, let 
R~ , R + , b, b\ be as in Lemma ^. Let a' p be the polynomial 
consisting of all the factors of a p whose spread is contained 
in W. Then 



d 



n 



(5) 



is a denominator bound of (1) with respect to W . 

Proof. Let y — p/q G K(n) be a solution of (1) and 
let u be an irreducible factor of q with multiplicity m and 
Spread(u) C W. We have to show u m \ d. Lemma 2 applied 
to p and v' implies that there is some i G ZT with u | q and 
u := N'u I a p . By the choice of s we have Disp 1 (u', u) < s. 

Lemma 4 implies the representation 



N p y 



Because of vi > 1, every i G R + differs from p in the first 
coordinate by more than s. This implies that N p u and 
hence that N p u m cannot appear in the denominator of N'y 
for any i G R + . But it does appear in the denominator 
of N p y, so it must appear as well in the denominator of 
the right hand side. The only remaining possibility is thus 
^ P " m I riiefl- N'^p, and hence 



n 

ie-R" 



JV 1 



Because of Spread(u') = Spread(u) C W, it follows that 
u m I d. □ 

The following figure illustrates the situation. The vector v 
is orthogonal to H but not necessarily to W , while the vector 
v' is orthogonal to W but not necessarily to H. Relation (4) 
separates p from the points in R + which are all below the 
plane H + sv. The points in R~ are all between H and 
H + sv. 




3.4 A Denominator Bounding Algorithm 

We now combine Theorems 2 and 3 to an algorithm for 
computing a denominator bound with respect to an arbi- 
trary given W in situations where these theorems are appli- 
cable. 

Definition 3. Let p, p' be corner points of S and W 
some submodule of TIT . 

1. The point p is called useless for W if there is an edge 
(p, s) in the convex hull of S \ {p} with p-s£ W. 

2. The pair (p,p') is called opposite if there is a vector 
v such that (s — p) ■ v > and (p' — s) ■ v > for all 
s G S. Such a v is called a witness vector for the pair 
(P,P')- 

3. The pair (p, p') is called useful for W if it is opposite 
and neither p nor p' is useless for W . 

The definition of a useful pair is made in such a way that 
when a change of variables as described in Section 3.1 is ap- 
plied which maps a witness vector of the pair to the first axis, 
then the sets A and B from Theorem 2 are such that p G A, 
p' G B (because of the oppositeness), no two elements of 
A agree in the first r — dimW coordinates (because p is 
not useless), and the same is true for B (because p' is not 
useless) . 

Whether a pair (p, p') G S 2 is useful or not can be found 
out by making an ansatz for the coefficients of a witness 
vector and solving the system of linear inequalities from the 
definition. The pair is useful if and only if this system is 
solvable, and in this case, any solution gives rise to a witness 
vector. 

If for a given submodule W we have found a useful pair, 
then we can compute a denominator bound with respect to 
W by the following algorithm. 

Algorithm 1. Input: An equation of the form (1), a 
submodule W of TIT , a useful pair (p, p') of S for W . Out- 
put: A denominator bound for (1) with respect to W. 

1 Sett := r - dim W. 

2 Choose vi, V2, . . . , vt G TIT such that vi is a witness 

3 vector for (p, p') and TIT = V ®W where 

4 V is the module generated by these vectors. 

5 Perform a change of variables as described in 

6 Section 3.1 such that Vi becomes the ith unit 

7 vector in WT , W becomes {0}* x % r ' t . 

8 Determine A, B as in Theorem 2. 

9 Compute s G N U { — oo} as defined in Theorem 2. 



10 Choose an inner vector v G IR' for p. 

11 Compute R~ as defined in Lemma 4- 

12 Compute d as defined in Theorem 3. 

13 Apply the inverse change of variables to d, getting d' . 

14 Return d' . 

The following variations can be applied for further im- 
provements: 

1. If the dimension of V is larger than 1, there might 
be different choices of witness vectors. Choosing dif- 
ferent versions in line 2 might lead to different de- 
nominator bounds of W, say, di, . . . ,dk. Then taking 
d :— gcd(di, . . . , dk) leads to a sharper denominator 
bound for (1) w.r.t. W . 

2. Choosing different inner vectors in line 10 might lead 
to different sets R~ to write (4) and hence gives rise to 
different denominator bounds in (5). Taking the gcd of 
these denominator bounds produces a refined version. 

We remark that the coefficients a s with s G S are of- 
ten available in factorized form. Then also the denomina- 
tor bounds are obtained in factorized form, and the gcd- 
computations reduce to comparisons of these factors and 
bookkeeping of their multiplicities. 

4. A COMBINED DENOMINATOR BOUND 

As mentioned earlier, when setting W = {0}, one is able 
to derive an aperiodic denominator bound for equation (1). 
In this particular case, for each corner point p there is an 
other corner point p' such that (p, p') is useful for W. Hence 
applying Algorithm 1 for any useful pair leads to an aperi- 
odic denominator bound. In particular, running through all 
corner points and taking the gcd for all these candidates 
leads to a rather sharp aperiodic denominator bound for 
equation (1) which coincides with the output given in our 
previous investigation [11]. 

In the other extreme, when setting W = TTT , a denomina- 
tor bound for (1) w.r.t. W would lead to a complete denom- 
inator bound for equation (1). However, in this case, we will 
fail to find a useful pair (p, p'), and our Algorithm 1 is not 
applicable. 

Our goal is to find a simultaneous denominator bound 
with respect to all W to which Algorithm 1 is applicable, 
i.e., for all W from the set 

U := {W submodule of TL r | 3 (p,p') useful for W}. 

In general, this is an infinite set. But we can make use of 
the observations made after Example 2. Using Lemma 2, it 
turns out that instead of looping through all these infinitely 
many modules W, it is sufficient to consider those W which 
appear as spread of some factor in the coefficient of a p . 
This argument even works for all W in the larger set 

O ■— {W submodule of IT 

3 (p, p') opposite with p not useless for W}, 

but since the W G O \ U do not satisfy the conditions of 
Theorem 2, we can only obtain partial information about 
their denominator bounds. 

We propose the following algorithm. 



Algorithm 2. Input: An equation of the form (1) . Out- 
put: A finite set of irreducible polynomials P = {pi, . . . ,Pk}, 
and a nonzero d G K[n] such that for every solution y = - G 
K(n) of (1) and every irreducible factor u of q with multi- 
plicity m exactly one of the following holds: 

1. Spread(u) G U and u m \ d, 

2. Spread(u) G O \ U and 3 s G V ,p G P : N s u = p, 

3. Spread(u) £ O. 

1 d := 1 

2 P:={} 

3 C := {p G S : p is a corner point of S} 

4 for all q G C do 

5 forall it | a q irreducible do 

6 W := Spread(u) 

7 if W G 17 then 

8 Compute a denominator bound do w.r.t. W 

9 using an arbitrary useful pair for W. 

10 d := lcm(d, do) 

11 else if W G O then 

12 P~Pu{u} 

13 return (P, d) 

Theorem 4. 27ie polynomial d computed by Algorithm 2 
is a denominator bound with respect to any finite union of 
modules in U . 

Proof. Let W be in U and (p, p') be a useful pair with 
respect to W . Let y = p/q be a solution of (1) and u be an 
irreducible factor of q with multiplicity m and Spread(u) C 
W. We have to show that u m \ d. 

Since p is not useless, Lemma 2 implies that there is some 
i G ZT with N'u | o p . This factor is going to be investigated 
in some iteration of the loop starting in line 5. The polyno- 
mial do computed in this iteration is a denominator bound 
with respect to Spread(A rl u) = Spread(w) and hence a for- 
teriori a denominator bound with respect to W. It follows 
that u m | do d. 

This proves the theorem when W itself is in U . If W is 
only a finite union of elements of U, the theorem follows 
from here by Lemma 1. □ 

For W G O \ U, we can still apply Lemma 2 but Theo- 
rem 2 is no longer applicable. This prevents us from com- 
puting precise denominator bounds with respect to these W. 
However, using the set P — {pi, . . . ,pt} returned by the al- 
gorithm we can at least say that for every denominator q of 
a solution y = p/q of (1) there exist m G N and a finite set 
S' C TL r such that 

d Yl N s p m (6) 

p£P 

ses' 

is a multiple of every divisor of q whose spread is contained 
in some finite union of modules in O. Appropriate choices 
S' and m can be found for instance by making an ansatz. 
Note also that the set P is usually smaller than the set of 
all periodic factors that occur in the coefficients a s of (1). 
This phenomenon was demonstrated already in the second 
part of Example 2. 

Summarizing, some part of the denominator is out of 
reach, namely all those parts of the denominator w.r.t. the 
modules from 

{W submodule of % r 

V (p, p') opposite for W with p and p' useless for W}, 



some part of the denominator can be given up to possible 
shifts and multiplicities, and a big part of our denominator 
bound can be given explicitly by d. 

The following improvements can be utilized. 

1. As preprocessing step, one should compute an aperi- 
odic denominator bound for the equation (1) as de- 
scribed above. What remains is to recover the peri- 
odic factors. As a consequence, one can neglect all 
irreducible factors u which are aperiodic and one can 
apply Theorem 3 where all aperiodic factors are re- 
moved from the polynomials a' p . 

2. Choosing different useful pairs for a module W in line 9 
might lead to different choices of denominator bounds, 
and taking their gcd gives rise to sharper denominator 
bounds of (1) w.r.t. W . 

5. DISCUSSION 

Typically the set O will contain all the submodules of Z r . 
Only when the convex hull of S happens to have two parallel 
edges on opposite sides, as is the case in Example 2, then 
modules W parallel to this edge do not belong to O. The 
set U will never contain all the submodules of 7L r . Precisely 
those modules W which are parallel to an edge of the convex 
hull of S do not belong to U. Since the convex hull of S 
contains only finitely many edges, U will in some sense still 
contain almost all the submodules of 7L r . 

Depending on the origin of the equation, it may be that 
there is some freedom in the structure set S. For example, 
by multivariate guessing [10] or by creative telescoping [16, 
9, 14] one can systematically search for equations with a 
prescribed structure set. In such situations, one can try 
to search for an equation with a structure set for which U 
and O cover as many spaces as possible. 

If two equations with different structure sets are available, 
it may be possible to combine the two denominator bounds 
obtained by Algorithm 2 to a denominator bound with re- 
spect to the full space 7L T . 

Example 4. Consider the following system of equations: 

-(k + n + l)(2fc + 3n + l)y(n, k) 

+ (k + n + 4)(2/c + 3n + 3)y(n, k + 1) 

- (k + n + 2) (2k + 3n + A)y(n + 1, fe) 

+ (k + n + 5)(2fc + 3n + 6)y(n + 1, k + 1) = 0, 
(n 2 + n + l)(2k + 3n + 3)y(n, k + 1) 

- (n +5n + 7) (2k + 3n + 4)y(n + 1, fc) 

- (n +3n + 3)(2k + 3n + 8)y(n + 1, k + 2) 

+ (n +7n+ 13)(2fc + 3n + 9)y(n + 2,k + 1). 

Algorithm 2 applied to the first equation returns 

d = (n + k + l)(n + k + 2)(n + k + 3)(3n + 2k + 1) 

as a denominator bound with respect to any W except (g)Z5 
and (jjz. Applied to the second equation, it returns 

d=(n 2 + n + l)((n + l) 2 + (n + 1) + l)(3n + 2k + 1) 

as a denominator bound with respect to any W except (jJZS 
and ^jJZS. The least common multiple of the two outputs is 
a simultaneous denominator bound with respect to any W . 



Indeed, the system has the solution 

i 

(n + fc + l)(n + fc + 2)(n + fc + 3)(n :4 +n+l)((n + l) 2 + (n+l) + l)(3n + 2fe + l) ' 

There is no hope for an algorithm which computes for any 
given single equation a denominator bound with respect to 
the full space TIT . This is because there are quations whose 
solution space contains rational functions with no finite com- 
mon denominator. For instance, for every univariate polyno- 
mial p, we have that l/p(n + k) is a solution of the equation 

y(n + l,k)- y{n,k + 1) = 0. 

It would be interesting to characterize under which circum- 
stances this happens, and to have an algorithm which finds 
a denominator bound with respect to TIT in all other cases. 
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